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A new version of the separable kernel of the nucleon-nucleon interaction in the Bethe-Salpeter approach 
is presented. The phase shifts are fitted to recent experimental data for singlet and uncoupled triplet partial 
waves of the neutron-proton scattering with total angular momenta J = 0, 1. The results are compared with 
other model calculations. 
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1. INTRODUCTION 

The relativistic description of the interaction be- 
tween two nucleons can be carried out using the Bethe- 
Salpeter (BS) equation pj. The BS formalism allows 
a covariant description of various reactions including 
electromagnetic ones. The BS approach with the sep- 
arable kernel of interaction (as a relativistic analogue 
of the Lippmann-Schwinger equation with the separa- 
ble potential [2l[3]) demonstrates a good description of 
electromagnetic processes with the deuteron (such as 
,the elastic lepton-deuteron scattering, deuteron electro- 
and photodisintegration) and light nuclei (deep inelastic 
.scattering of leptons), see the review [i] and references 
therein. For instance, the BS formalism facilitates anal- 
ysis of the role of P waves (negative energy partial- wave 
components of the BS amplitude) in the electromagnetic 
properties of the deuteron and its comparison with the 
nonrelativistic treatment [5]. Furthermore, the covari- 
ant BS approach makes it possible to analyze off-mass- 
shell effects and contributions of the relativistic two- 
body currents. 

The consistency of the consideration of the deuteron 
breakup reactions demands the final state interaction 
(FSI) between the outgoing nucleons to be taken into 
account. To introduce the FSI of the final np pair 
the BS equation for continuous state should be solved. 
These calculations were performed in the paper [B] for 
the separable kernel of nucleon-nucleon (NN) inter- 
action with a relativistic generalization of Yamaguchi 
functions. The obtained kernel allows one to describe 
the NN phase shifts up to laboratory kinetic energy 
^Lab^ 0.6-0.7GeV (for example, the nonrelativistic 
Bonn potential can be applied up to TLab^ 0.35 GeV). 
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However, when the amplitudes of reactions in a high en- 
ergies region are calculated, this type of functions leads 
to the presence of nonintegrable singularities. Therefore 
there are many approaches called quasipotential where 
the two-body Green's function in the BS equation is 
replaced by a suitable function which allows to reduce 
the 4-dimensional integral to the 3-dimensional one, like 
in[7j. 



We construct a relativistic generalization of the mod- 
ified Yamaguchi form factors. For scalar nucleons such 
functions were introduced in [5] . Functions of this type 
do not contain poles on the real axis of the relative en- 
ergy. We use them for the description of the proton-neu- 
tron scattering in the wide region of TLab(till 3 GeV). It 
is necessary when we consider reactions with light nuclei 
with a high momentum transfer. In particular, the cal- 
culations of the (i(e, e'p)n reaction are supposed to be 
performed for the relativistic energies using the separa- 
ble functions elaborated in the paper (see, for example, 
different kinematical conditions in the Table from [S]). 



In this work the one-rank kernel of the NN inter- 
action for singlet and uncoupled triplet partial states of 
the neutron-proton (np) scattering with total angular 
momenta J = 0, 1 are presented. The paper is orga- 
nized as follows: in Sec. 2 we set out basic formulae of 
the formalism, in Sec. 3 the functions under considera- 
tion are described, in Sec. 4 we show the scheme of the 
calculations and results. Sec. 5 is devoted to the analy- 
sis of the results and the comparison with calculations 
within other models, and in Sec. 6 we summarize the 
approach. 
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2. BASIC FORMULAE 

Let us consider the partial-wave decomposed BS 
equation for the nucleon-nucleon T matrix in the rest 
frame of two particles (the square of the total mo- 
mentum s = (pi + ^2)^ and the relative momentum 
p = {pi —p2)/2 [p' — (p[ —^2)72] are defined through 
the initial pi, p2 and faialpi, p'^ nucleons 4-momenta): 

THip'o,p';PQ,p;s) = Vu{p'„,p';po,p;s) 

-I-— / dko / k'^dkVii{pQ,p';ko,k;s) 

xS{ko,k;s)Tu{ko,k;po,p;s). (1) 

Here Tu is the partial-wave decomposed T matrix (I 
enumerates ^Lj states for simplicity, S is the spin of 
two particles, L is an orbital momentum, and p spin de- 
fines a sign of an energy), Vu is the interaction kernel, 
Ek = \/W^+rn? , ?7i is a nucleon mass. We consider un- 
coupled partial states, and therefore the T matrix does 
not have transitions to different partial states. The two- 
particle propagator is 

S{ko,k- s) = 1/ ((Vi/2 -Ek+ ief - kl) . (2) 

Using the separable anzatz for the one-rank kernel 
of the NN interaction 

Vii{pq,p';pq,p;s) = Xi{s)g^^\pQ,p')g^^\po,p), (3) 

we obtain the solution of the equation ^ for the T 
matrix in the similar separable form: 

Tnip'o,p';po,p;s) = Tiis)g^'\p'o,p')g^'Kpo,p), (4) 

where 

Ti{s)^l/{Xi{s)-^+hi{s)) (5) 

with 



hiis)=-^ I dko I k^dk 



{g^Hko,k)r 



(Vi/2 -Ek+ i£)2 - k^ 



.(6) 



and A; is a model parametric function. From the nor- 
malization condition for the T matrix in the on-mass- 
shell form: 



Tu{s) = Tu{0,p;0,p;s) 



167r 



JSi 



y/sy/s — 4™^ 



e^'^ sin Si, (7) 



where p = •\/s/4 



•^/mTLab/S is the on-mass- 



shell momentum, the phase shift Si for the one-rank ker- 
nel of interaction is defined as 

cot Si{s) = ———— = — . (8) 



We obtain the low-energy scattering parameters ex- 
panding the following expression for S waves into series 
of p-terms [lOj: 

pcot^,(.s) = -1/ai + ^p'roi +0(p3). (9) 

Here the scattering length ai and the effective range roi 
are introduced. Since it is clear which partial states are 
discussed, hereinafter we omit the index I for simplicity. 



3. ONE- RANK KERNEL 

We analyze two relativistic generalizations of the Ya- 
maguchi form factors: modified Yamaguchi (MY) func- 
tions and modified extended Yamaguchi (MEY) func- 
tions. 

31 Modified Yamaguchi functions 

For the description of the chosen partial states we 
use the following covariant expressions: 



9^'Hpo,p) 



.9'^1(P0,P) 



(Pc 

(p2_p2_^2)2 + 

V-Po+P"^ 



(10) 



(11) 



ipl-p^~l3^y+a^' 

The functions are numerated by angular momenta L = 
0([S']), 1([P]). The numerator in g^^^ is introduced to 
compensate an additional dimension in the denomina- 
tor to provide the total dimension as GeV~^. This form 
was chosen because at pd = /3^ , a = we get the func- 
tion gl^l in the standard Yamaguchi form p. We prefer 
not to consider the case with the square root in the de- 
nominator as it is used in [8] for 5* waves to avoid the 
calculations with cuts on the real axis in pq complex 
plane. 

32 Modified extended Yamaguciii functions 

To extend the form of functions with increasing num- 
ber of parameters we introduce the following g func- 
tions: 



JS] 



iPo,p) 



g^'^Kpo.p) 



(PcI-pI+P^) 
ipl-p'-Pir + ai 

Ci2{pI-P^){Pc2-pI+p'? 
{{pl-p^-Pl?+alY ' 

\/-pI + P^ 



(12) 



aj 



[pl-p '-0if 

, Cxi Vi-Po+P^)^iPc3 -pI+ p") 



Im T/i(s) 



lmhi{s) 



{{pl-p2-pl)2+al) 



vi\2 



(13) 
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The denominators with pci , Pc2 , Pc3 are intended for di- 
mension compensation as pd for g'- ' in the previous 
case (fTOl). 



4. CALCULATIONS AND RESULTS 

Using the np scattering data we analyze the parame- 
ters of the separable kernel distinguishing three different 
cases: 

1. There are no sign change in phase shifts or bound 
state {^Pi, ^Pi partial states). In this case 



Xi{s) — const. 



(14) 



for the phase shifts till the maximal laboratory ki- 
netic energy rLab< Thab^^ deuteron energy and low- 
energy scattering parameters. Experimental data 
for phase shifts are taken from the SAID program 



(http://gwdac.phys.gwu.edu) and for the rest quanti- 



ties - from [TT]. For all states except '^Pq the Tj 



Lab 



is 



chosen to be '^ 1 GeV both for MY and MEY cases and 

for ^Po wave 3 GeV for MEY case. 

Now we find the introduced parameters of the kernel: 

1. For P waves the minimization procedure for the 
function 

n 

X' - E(<5'='=P(s.O - S{s,)f/{A6'^'^^{s,jf (19) 



This is sufficient for the most of the higher partial 
waves. 

2. One sign change and no bound state (^S'o, ^Po 
partial states). In this case the energy-dependent 
expression for A; is used (see 8 and references 
therein) : 



A/(s) = A(so - s). 



(15) 



Here the parameter sq is introduced to reproduce 
the sign change in the phase shifts at the position 
of the experimental value for the kinetic energy 
TLab where they are equal to zero. It is added to 
the other parameters of the kernel. 

3. One sign change and one bound state {^Si state). 
In this case in addition to the zero in the phase 
shifts the T matrix has a pole at the mass of the 
bound state M^: 



det|r-i(s = Mj2)| =0 



To achieve this we set 



Ms) = A- -2, 



s — mf, 



(16) 



(17) 



The parameter toq is chosen to satisfy the condi- 
tion (see eq. p6)) and pT| l 



So - M 



-1 ^^^ 



'0 



Reh{M^) =0, 



(18) 



where Mf, = {2m — Eb) defines the mass of the bound 
state and Eh is its energy. 

The calculation of the parameters is performed by 
using the equations dH), ([9|) and expressions given 
in previous section to reproduce experimental values 



is used. Here the number of experimental points 
n depends on the T"^if value. 

2. For S waves the values of the scattering length a 
are also included into the minimization procedure 

n 
;^2^^(5-P(.,)-^(s,))V(A<5°'=P(..))' 



-(a'^"P-a)V(Aa""P)'. 



(20) 



The effective range tq is calculated through the 
obtained parameters and compared with the ex- 
perimental value Tq'^^. The value of the binding 
energy E^i for the case ^Si is taken into account 
by introduction of the parameter m^ in eq. (|17p . 

The calculated parameters of the one-rank kernel 
with MY and MEY functions and T^^ are listed in 
Tables 1 (here the values of Sq are presented, too) and 2. 
In Table 3 the calculated low-energy scattering param- 
eters for S waves are compared with their experimental 
values. 

In Figs[T]-[5] the results of the phase shifts calcula- 
tions are compared with experimental data and two al- 
ternative descriptions by CD-Bonn [12^ and SPOT [13,. 



5. DISCUSSION 

The calculated phase shifts are presented in Figsd]- 
[5] and parameters of the separable kernel are given in 
Tables 1-2 (MY case is denoted by the dashed line, 
MEY one - by the solid line). The calculations with 
the nonrelativistic CD-Bonn potential (dotted line) and 
the solution SP07 [13] (dash-dotted fine) are included 
for comparison. 

We can see in FiglT] that all the calculations give 
a reasonable agreement with experimental data up to 
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energy TLab~ l.lGeV but at larger energies their be- 
havior becomes different drastically. Note that the CD- 
Bonn and SP07 results for phase shifts even change 
sign at 1.5 < TLab < 2.5 GeV. So we can stress that 
it is very desirable to obtain more exact determination 
of experimental data for the ^Pi channel in this en- 
ergy range. Now keeping in mind that the discrepan- 
cies between CD-Bonn, MEY and SP07 calculations at 
TLab < 1-2 GeV are not so large we can trust our calcu- 
lations in this range of energies. 

In Fig[2]we can see that CD-Bonn and MY calcula- 
tions demonstrate an opposite behavior for '^P^^ phase 
shifts at TLab > 1 GeV. The MEY and SP07 calculations 
give a very good agreement with experimental data in a 
wide range of energies TLab < 3 GcV. Thus, our model 
for the kernel of the NN interaction in the ^Pq channel 
is acceptable to be used in various relativistic calcula- 
tions of reactions. 

Calculations for the "^P^ channel (see Fig 13]) show us 
that MY, CD-Bonn, and MEY have similar behavior in 
a wide energy range, but can explain experimental data 
up to TLab =0.55, 0.6 and 1.2 GeV, respectively. The 
SP07 calculations give a suitable agreement with exper- 
imental data up to 3 GeV. Thus to explain the data for 
the "^P^ channel in wider energy range in our approach 
we need to increase the rank of the separable kernel of 
the NN interaction. 

FigH] demonstrates a close agreement with exper- 
imental data for CD-Bonn up to the energy TLab < 
0.6 GeV, MY and MEY up to TLab < 1.2 GeV in the 
^Sq channel (note that these calculations give very sim- 
ilar behavior of phase shifts in the wide energy range up 
to 3 GeV). We can conclude that in this case as well as 
for ^P^ we need to increase the rank of the separable 
kernel of the NN interaction. The SP07 calculations 
describe experimental data for the ^Sq channel up to 
3 GeV. 

FigOfor the "^S^ wave shows us that CD-Bonn cal- 
culations agree with experimental data up to 0.6 GeV. 
MY, MEY, and SP07 have similar behavior in a wide en- 
ergy range and explain known experimental data up to 
TLab =1.1 GeV. It should be noted that calculations of 
low-energy parameters with MY and MEY form factors 
give us a reasonable agreement with their experimental 
values (see Table 1). 



MEY functions give a better description of the scatter- 
ing data than the MY ones. Using the suggested MEY 
form factors the phase shifts are described in a whole 
range of measured energies for the following partial 
states: iF+, ^S+ (TLab < 1.2 GeV) and ^P+ (TLab < 
3 GeV). The phase shifts for the ^S^ and 3p+ partial 
states can be described in our approach up to TLab < 
1.2 GeV. To improve an agreement with experimental 
data up to TLab < 3 GeV it is necessary to increase the 
rank of the separable kernel of the NN interaction. It 
is planned to be done in the near future. 
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6. CONCLUSION 



We present the new parametrizations of the separa- 
ble kernel of the NN interaction which are adopted for 
calculations at large energies. As it was expected the 
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Parameters of the one-rank kernel with modified Yamaguchi functions. 





lc+ 
Do 


35+ 




ip+ 


3p+ 


3p+ 


7^ in ax 
-'Lab 


0.999 GeV 


1.1 GeV 




1.1 GeV 


0.999 GeV 


0.999 GeV 


A (GcV2) 


-1.0375 


-10.004 


A (GeV^) 


8.11284 


-211.74 


11.962 


a (GeV) 


1.3312 


-11.287 


13 (GeV) 


0.56487 


0.53808 


0.58466 


/3 (GeV) 


0.1 


1.3540 


a (GeV) 


0.2 


0.85720 


0.2 


Pel (GeV2) 


19.229 


0.49696 










so (GeV2) 


4.0279 


4.2020 


3.8682 



Parameters of the one-rank kernel with modified extended Yamaguchi functions. 





lc+ 

Dq 


35+ 




ip+ 


3p+ 


3p+ 


-'Lab 


0.999 GeV 


1.1 GeV 




1.1 GeV 


2.83 GeV 


0.999 GeV 


A (GeV2) 


-0.84694 


-9.1434 


A (GeV^) 


0.054821 


-192.60 


0.030177 


Cu (GeVO) 


-8.6404 


3.3641 


C12 (GeVO) 


4.67839 


4.9604 


-0.39145 


/3i (GeV) 


0.13400 


2.6609 


Pi (GeV) 


0.21829 


0.71329 


0.19951 


(32 (GeV) 


0.92330 


0.12501 


(32 (GeV) 


0.20829 


0.70329 


0.66216 


ai (GeV) 


1.4553 


2.0808 


ai (GeV) 


0.20500 


0.70958 


0.20000 


a2 (GeV) 


0.93129 


1.6340 


a2 (GeV) 


0.52835 


1.4574 


0.32015 


Pel (GeV2) 


24.261 


189.85 


Pes (GeV2) 


0.36868 


-1.7576 


-117.75 


Pe2 (GeV2) 


-2.1076 


3.6688 











The low-energy scattering parameters and binding energy for the singlet (s) and triplet (i) S waves. 



a,(Fm) ro,(Fin) at{Fm) rot(Fm) Ed{McY) 



MY 


-23.80 


2.4 


5.43 


1.8 


2.2246 


MEY 


-23.74 


2.51 


5.39 


1.73 


2.2246 


Experiment 


-23.748(10) 


2.75(5) 


5.424(4) 


1.759(5) 


2.224644(46) 
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T ,GeV 

Lab 



Phase shifts for the ^P^ wave. Dashed hne corresponds to our parametrization with MY functions, ea. (|lip : the 
solid line illustrates the extended case, ea. (|13p : the dotted and dash-dotted lines describe the CD-Bonn and SP07 
calculations, correspondingly. 
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CO 

to 




-100 - 



SP07 



T ,GeV 

Lab 



The same as in Fig.l for the Pq state. 
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-60 - 



-70 




T ,GeV 

Lab 



The same as in Fig[T]for the P^ state 
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T ,GeV 

Lab 



The same as in Fig[T]for the singlet partial state ^Sq, but the description by MY functions is defined by eq. (fT(7)l 
and extended one (MEY) - by ea. (fT2|l . 
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150 
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50 - 
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CO 



-50 - 



-100 - 



-150 




T ,GeV 

Lab 



The same as in Fig|4]for the triplet partial state S^ 



3q+ 



